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Abstract 

We revisit the general topic of thermodynamical stability and critical phenomena 
in black hole physics, analyzing in detail the phase diagram of the five dimensional 
rotating black hole and the black rings discovered by Emparan and Reall. First we 
address the issue of microcanonical stability of these spacetimes and its relation to 
thermodynamics by using the so-called Poincare (or "turning point") method, which 
we review in detail. We are able to prove that one of the black ring branches is 
always locally unstable, showing that there is a change of stability at the point where 
the two black ring branches meet. Next we study divergence of fluctuations, the 
geometry of the thermodynamic state space (Ruppeiner geometry) and compute the 
appropriate critical exponents and verify the scaling laws familiar from RG theory 
in statistical mechanics. We find that, at extremality, the behaviour of the system is 
formally very similar to a second order phase transition. 
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1 Introduction 

Topology change of event horizons constitutes a fundamental issue in gravity. From a 
broad perspective, one of its most important implications comes from its relation to the 
problem of resolution of singularities. Instabilities leading to the destruction of an event 
horizon can in principle produce naked singularities, hence violating the Cosmic Censorship 
Conjecture. The latter constitutes one of the major unsolved issues of the classical theory, 
and a complete understanding of this problem must imply deep consequences at the quan- 
tum level. Somehow, the Cosmic Censorship can be thought as the classical mechanism 
for resolution of singularities, which ultimately are expected to be cured by any correct 
quantum completion of gravity. 

In a String Theory context, singularity resolution involving a (manifold) topology 
change was first emphasized in P, where the inclusion of nonperturbative string states 
(namely, black hole states) was shown to provide a mechanism for the resolution of the 
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conifold singularity. Gravitational phase transitions also play an important role in the 
AdS /CFT correspondence, starting from the original work of Witten [2] on nucleation of 
black holes in Anti-de-Sitter space (the Hawking-Page phase transition jH]) and its rela- 
tion to deconfinement at finite temperature in gauge theories. This is of interest both for 
the study of the dual phenomena in the corresponding gauge theory (recent works include 
e.g. jl]) and, also, because it seems to provide the possibility of investigating gravitational 
singularities from the field theory point of view [S]. 

Interestingly, higher dimensional General Relativity provides indications of topology 
change of event horizons, the best known example being the Gregory-Laflamme instabili- 
ty jnj. Recently this phenomenon has been studied in the presence of compact dimensions 
by a number of authors [7j (for a recent review see 0). In general, the question of topol- 
ogy change of event horizons makes sense whenever different phases of black objects (black 
holes, black strings, black p-branes or others) can exist, a general feature of higher dimen- 
sional General Relativity and low-energy String Theory. 

Of course, and at least at classical level, the key issue to be investigated about possible 
transitions between different black phases concerns their dynamical stability, i.e. stability 
of the solution under field perturbations, or classical stability under Lorentzian time evo- 
lution. However, not much is known about the classical stability of solutions of GR above 
four dimensions. Apart from the well known instabilities of the Gregory-Laflamme type, 
very little is known about many other solutions of interest, like e.g. those with angular 
momentum. For those, some partial or indirect checks of dynamical stability have been 
put forward 9j but, in most cases, an exhaustive and rigorous proof of (in)stability is still 
lacking. 

A way to address the problem of stability of event horizons comes from black hole 
thermodynamics. In ordinary thermodynamics of extensive systems, local thermodynamical 
stability (defined as equivalent to the condition that the Hessian of the entropy has no 
positive eigenvalues) is related to the dynamical stability of the system. However, such a 
parallel cannot be inferred for black holes, which constitute non-extensive, non-additive 
thermodynamical systems. The reasons for this depart from standard thermodynamics 
can be phrased in different ways: the long-range behaviour of gravitational forces, the well 
known scaling of the BH entropy with the area (rather than the volume), or the fact that 
BHs constitute elementary entities that cannot be subdivided into separate systems, not 
even in any idealized manner. 

Critical thermodynamic behaviour of black holes was first studied in and, since 
then, by a number of authors (see e.g. ^T] and references therein). However, the interpre- 
tation of the divergences in the specific heats and generalized susceptibilities (in particular, 
their interpretation as a true phase transition and their relation to the stability of the sys- 
tem) is not clear and has been subject of debate [121 HSl UH Uil UHl UIl CHI- In a String 
Theory context some partial attempts to relate dynamical stability and local thermody- 
namical stability have been put forward as well but currently the present status of 
this subject is not clear jSHI (see also the third reference listed in |3]). We will comment 
on the relation between the Gubser-Mitra conjecture and the "Poincare method" that we 
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will use in this paper in the Conclusions. 



In this paper we will explore these issues in the context of five dimensional gravity. This 
scenario seems interesting to us due to the discovery made by Emparan and Reall of a new 
BH phase: an asymptotically fiat, rotating black hole solution with horizon topology x 5"^ 
and carrying angular momentum along the — what they called a "black ring" [21] . Such 
black ring solutions exist in a region in parameter space (i.e. mass and angular momentum) 
which overlaps with that of the five dimensional spherical rotating black hole with rotation 
in just one plane studied in |221- This fact provides the first known example of BH non- 
uniqueness in asymptotycally fiat space, and it raises questions on stability and on the 
possibility of transitions between the different phases. 

Let us summarize the main properties of these spacetimes. The five dimensional rota- 
ting black hole |22| has an upper "Kerr bound" in its angular momentum per unit mass, 
the extremal configuration that saturates the bound being a naked singularity. In what 
follows, we will be using the dimensionless quantity 

2 

X ^ 

(where J and M are respectively the total angular momentum and mass) as the "control 
parameter"^ of the problem — analogous, say, to the temperature in a liquid-gas system 
when studied in the canonical or grand canonical ensembles. The same parameter has been 
already used in investigations of black rings in [22] ■ It is normalized in such a way that 
the five dimensional rotating black hole solution exists in the range 

Black Hole : < a; < 1 , 

X = 1 being the extremal (singular) limit. On the other hand, one can have two different 
black ring spacetimes with horizon topology x S'^ when the total angular momentum 
per unit mass exceeds a certain minimal value Xmin given by: 




^m\n \ ~ 0.919 



One of them, that we will call the "large" black ring, can have angular momentum per 
unit mass unbounded from above. The other one, that we call "small" black ring, cannot^. 
One has: 

Large Black Ring : 

Small Black Ring : Xmin < a; < 1 . 

At Xmin both solutions coincide. For a comparison of the respective entropies see Fig. ^ 
One can see that near extremality the large black ring is entropically favoured (its en- 
tropy equals that of the rotating BH at x = 2a/2/3 ~ 0.943). This fact was taken in 

proper order parameter will be defined in Section |31 
^The labels "large" and "small" for the two different black rings have already been introduced in |23|. 
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Ref. 21] as a first indication of a pliase transition from tlie black liole to a black ring 
when the former rotates fast enough, i.e. as we vary x. Even if such a topology changing 
process is not believed to be possible classically plj, one might expect to find some evi- 
dence of it by exploring issues like stability and the thermodynamics of these spacetimes. 
The study of the general properties of this phase diagram is the main purpose of this paper. 

There are two main concrete questions that we will address here. The first one is that 
of the classical stability of these solutions and its relation to thermodynamics. We will deal 
with this issue in Section 13 and our main tool will be the so-called "Poincare" or "turning 
point" method. 

The second question is that of to which extent we are allowed to call all these configura- 
tions "phases" in any thermodynamical or statistical mechanical sense. This is a nontrivial 
question which has to do with the properties of the classical "configuration space" (i.e. the 
space of metrics): namely, if all these solutions correspond or not to extrema of some com- 
mon underlying thermodynamic potential. If so, the general theory of critical phenomena 
tells us that, near a critical point, the generalized susceptibilities and fluctuations have to 
diverge in a very specific way. We will address these issues in Sections |3] and El Our main 
tools will be the study of the geometry of the thermodynamic space state and the general 
theory of second order phase transitions. 

Summary of Results 

After reviewing in Section El the main properties of rotating black holes and black rings, 
we focus in Section (21 on the issue of stability. We revise the usual arguments relating 
dynamical and local thermodynamical stability, and recall why these cannot be applied to 
non-extensive systems such as BHs. Instead, we will use the so-called "Poincare method 
of stability" j2S| as it has been already applied to BH physics in [THj ITT)! IT7| EE] ■ This 
technique is of remarkable simple applicability, and it essentially consists on plotting the 
appropriate phase diagrams. We will review it in detail and comment on the suitability 
of this method for the study of BHs. Contrary to the standard analysis based on the 
sign of the specific heat, the Poincare method does not predict any instability in e.g. the 
Schwarzschild or Kerr spacetimes. We will use it to study the BH/BR system in the 
microcanonical ensemble, finding that the only point at which a change of stability can 
occur is at Xmin- One of our main results is in fact that the small black ring is not only 
globally, but also locally unstable. The large black ring turns out to be a "more stable" 
phase in a sense that we will explain. On the other hand, we do not find any change 
of stability along the whole BH branch. In particular, we will see that nothing special 
happens at the point in the BH branch where the specific heat changes sign (that is, the 
five dimensional analogue of the "critical point" found by Davies in jTO] and discussed since 
then by many other authors). 

In Section |31 we discuss divergence of fluctuations and critical behaviour. We will see 
that quadratic fluctuations can be evaluated at x^in and x = 1, due to the fact that at 
those points the relevant degrees of freedom of the system are just a few and under control. 
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We will see that at both points fluctuations diverge, but that the origin of this divergent 
behaviour is very different at Xmin and at extremality. Next we will use the geometry of 
the thermodynamic state space as introduced by Ruppeiner [22], and explain how this 
technique should be interpreted and applied to the BH case. In particular, the divergences 
of the "thermodynamic curvature" are interpreted in Ruppeiner theory as indicators of 
critical phenomena. We will compute this quantity for the BH/BR system, finding that it 
diverges both at Xmin and a; = 1. 

In Section |S1 we compute the various critical exponents at the points where fluctuations 
diverge. We check the usual scaling laws familiar from Renormalization Group theory in 
statistical mechanics, finding that these are obeyed at extremality, but not at Xmin- We 
thus conclude that x = 1 might be identified, in a certain sense that we will explain, with 
the critical point of a second order phase transition. The phenomenon at Xmm seems, on 
the other hand, very different: there fluctuations diverge because, as a consequence of 
the change of stability, the thermodynamic potential becomes locally flat. We discuss our 
results and prospects for future work in Section IHl 



2 The Black Hole and Black Ring Spacetimes 
2.1 The d = b Rotating Black Hole 

In [2^ Myers and Perry found black hole solutions of higher dimensional GR with arbitrary 
angular momenta. In particular, the d = 5 BH solution with rotation in a single plane is 
found to be given by the line element: 

ds"^ = -dt^ + A{dt-a sin^ ed^f + sin^ 9 (r^ + a^)#2 + ^ c/r^ + dO"^ + cos^ Odcj)'^ , (2. 1) 

with < ^ < 27r, < < 27r, < 6* < 7r/2 and the functions A, \E' and p being: 



A^^-J^ p^^r^ + a\os'e. (2.2) 

+ cos^ o r"^ + a"^ — /j, 



The parameters /i and a set the physical ADM mass and angular momentum via the 
relations: 

M = ^/i, J = ^Ma. (2.3) 

This solution has a "ring singularity" at p = and a event horizon (of spherical topology) 
ai rjj = fi — , whose existence thus requires the bound 

n 32G , 

The extremal case (p = a^) is a naked singularity. 
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Thermodynamic Quantities 

Let us introduce the following dimensionless quantity (a suitably normalized angular mo- 
mentum "per unit mass"): 

_ ( 21ttV^^ J _ a 

In terms of it and the mass parameter /i, the area, temperature and angular velocity of 
the horizon are given by: 



^BH = 27rV'/'v/r^, T^^ = -^Vl^, fiBH = ^. (2.6) 

In particular, one can see that the area decreases when ^ at fixed mass (i.e. x — *■ 1), 
and it vanishes precisely in the extremal limit. This is because, as we increase the angular 
momentum of the BH, the event horizon gets "flattened" , becoming a (singular) disc of 
zero area in the extremal limit (for an explicit analysis of this geometry see e.g. 0). 



2.2 The Black Ring Spacetime 

Let us review now the solution found by Emparan and Reall in [21]^. This solution describes 
a five dimensional, asymptotically flat black hole whose horizon has x S'^ topology. The 
solution is regular on and outside the event horizon provided that it has angular momentum 
along the direction. Both the black ring and the rotating black hole solution are 
described by the metric: 



ds^ = y-idt + RV\i^il + y)diP] + 

Fiv) V / 



ny) 

K 



2 



{x-yy 



F{y) 2\ , N / dx^ , G{x 



F{x) G{y) d^' + -l!^dy'\+ F\y) -— + 



G{y) ^ y \G{x) F{x) ' 

where the functions introduced above are: 

F(0 = l-Ae, G(0 = (l-a(l-^0- (2 



(2.7) 



R and v are parameters whose appropriate combinations give the total ADM mass and an- 
gular momentum, while A has to take specific values in order to avoid conical singularities. 
In particular, there are two possibilities: 

1 , (black hole) , 

2^ (^.^ ^. ■ ^ ^^'^^ 

(black rmg) . 



1 + 



^Note however that here we will be using the conventions of |28| . 
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Once A is fixed one has therefore two parameters: R (with dimensions of length and setting 
the "radius" of the BH or the black ring) and v (dimensionless and varying from to 1). 
In the BH case = represents a non rotating black hole, while v = 1 yields the extremal 
case. In the BR case, and as explicitly shown below, as v takes values from to 1 we 
get two different black ring configurations, which coincide only ati^ = l/2. Atz/ = 
we have an infinitely thin and large black ring with infinite angular momentum per unit 
mass. The angular momentum per unit mass of this object takes its minimal value at 
V = 1/2. From = 1/2 to = 1 the solution describes a second black ring of different 
geometry, which approaches the singular BH spacetime in the extremal limit v ^ 1. At 
u = 1/2 both solutions coincide, yielding the black ring spacetime with minimum possible 
angular momentum per unit mass. These properties can be visualized e.g. from the entropy 
diagram of Fig. Q 

In j2HI the physical parameters of both the black ring and black hole solutions are found 
to be: 

,,3VA(1 + A)(A-z/)3/2 37ri?2A(A + l) 



A = STr'n-" ' — , M 



(l + z/)2(l -z/) ' AG u + l 



T = — J = v^v'-v^-r^v; (2.10) 

47ri? ^A(A - v) ' 2G (1 ■ • 



z/)2 



1 I V 



RV A(1 + A) ■ 



We recover the black hole or black ring values by using ()2.9p . Finally, let us recall that the 
physical parameters in ()2.10j) satisfy the Smarr relation given by 

M =- C^ + Qj] (2.11) 



2\4G 

for any value of A, which in particular means that both the black hole and the black ring 
spacetimes satisfy the same Smarr relation. 



Thermodynamic Quantities 

In order to compare entropies and eventually carry out a thermodynamic analysis, we will 
need a fundamental relation A = A{M, J) and the corresponding "equations of state" for 
the black ring spacetime too. Defining the parameters fi and x in terms of the physical 
mass and angular momentum as in ()2.3|) and ()2.5|) . one finds: 

^BR = v^vrV^/^v/Kr^, Tbr = — =^./^, i7BR = ^./^, (2.12) 

where u = z/( J^/M^) is given by any of the following two values: 

Us = \/8 X { 4= cos 9 + sin 9 ) — 1 , 

^ > { (2.13) 

i^L = v8 X ( 4= COS 9 — sin 9 ) — 1 , 
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the angle being 



1 / 32 

e = -arctanW — x2 - 1. (2.14) 
3 V 27 ^ ^ 

These exphcit expressions are obtained by inverting the parametrization from ()2.10|) and 
taking into account the conditions to avoid naked singularities^. The fact that the angular 
momentum at fixed mass is bounded from below at at the precise value Xmin = a/27/32 ^ 
0.919 [21^ is explicit in the formula above. 

In ()2.13|1 Us and ul correspond to the two possible black rings found in j2Tj: i^s cor- 
responds to what we will call "small" ring (i.e. the one whose angular momentum is also 
bounded from above at a; = 1), while z/^ corresponds to the "large" black ring (the one with 
unbounded angular momentum). As explained, in the limit x ^ oo the horizon of this 
large black ring takes the shape of an infinitely large (along the S^) and thin ring, hence 
the name. The fact that the angular momentum of the small BR is also bounded from 
above, while the one of the large BR remains unbounded, can be deduced by inspection 
of ()2.13p and recalling that u is constrained to take values in [0, 1) 



2.3 Entropy Diagrams 

The expressions ()2.(ij) and ()2.12j] allow for an explicit comparison of the respective entropies. 
We recover the plot already found in [21], which is shown in Fig. ^ 




Figure 1: Plot of A/2'K^fj,^/^ against x = ^JTJ-k ji^GJ jW^I"^ for the Myers-Perry BH (dashed hne) and 
the two black rings (see in])- At Xmin = a/27/32 ~ 0.919 black ring formation becomes possible. At 
X = 2^2/3 0.943 the entropy of the large BR exceeds that of the black hole. 



This figure is the first indication of the possibility of a phase transition when we vary 
the angular momentum at fixed mass. As we see, beyond x = 2a/2/3 ~ 0.942 the entropy 
of the large black ring exceeds the one of the black hole, which rapidly decreases to zero at 
X = 1. Emparan and Reall pointed out in their original work that this fact may be taken 
as an indication of a BH/BR phase transition. 

"^This fundamental relation has already been given in parametric form in '2'6'. 
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Motivated by this fact, we will study the stability properties and the behaviour of 
fluctuations of these spacetimes. However, our results will lead us to the conclusion that, 
from the point of view of stability, the relevant points in this phase diagram are Xmin and 
X = 1, and that nothing special seems to be happening at x = 2v^/3. We will comment 
on the implications of this result in the Conclusions. 



3 Dynamical vs. Thermo dynamical Stability 

To argue in favour or against a possible phase transition, the first issue we should address 
is that of the stability of the classical solutions. To face this problem, one may think that 
a way to circumvent the study of linear perturbations of the metric might come from ther- 
modynamics. In ordinary thermodynamics of extensive systems, local thermodynamical 
stability (defined as the Hessian of the entropy having no positive eigenvalues) is linked 
to the dynamical stability of the system: a positive mode in the Hessian means that at 
least some kind of small fluctuations within the system are entropically favoured, implying 
that the system is unstable against these. However, the simple argument leading to this 
conclusion [21] relies on the additivity of the entropy — a property which does not hold for 
black holes. The best known example of this failure is Schwarzschild black hole: a stable 
configuration with negative specific heat (i.e. a positive Hessian). 

In order to make clear the differences with the case of BH thermodynamics, let us 
first review the standard criteria of local thermodynamical stability. To fix ideas, consider 
an extensive thermodynamic system in the microcanonical ensemble with a fundamental 
relation S = S{M, J) (where M and J stand for the internal energy and some other 
extensive parameter). Divide it into two identical subsystems and consider a spontaneous 
fluctuation that transfers some amount of energy dM from one subsystem into the other. 
Then, if the system is to be stable, one requires the total entropy after this energy transfer 
to be smaller (or equal) than the previous total entropy: 

S{M + dM, J) + S{M - dM, J) < 2S{M, J) , (3.1) 

since, otherwise, energy transfers within the system will take place spontaneously, therefore 
meaning that the system would be unstable against redistributions of mass from homo- 
geneity (this would be the dynamical unstable mode). The differential form of the above 
condition is: 

ffMM-I^SO. (3,2) 
So ()3.2|) is a necessary condition for a system to be stable. Similar considerations lead to 

Hjj^^<0 (3.3) 
for transfers of J within the system. For "cooperative" transfers of both M and J we have 
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to require the Hessian H of S'(M, J) to not to be positive definite. Hence 



) 



2 



detH 



dM^ (9J2 



dMdJ 



> 0. 



(3.4) 



So conditions ()3.2|) - ()3.4j) have to be satisfied if the system is to be stable. Only two of 
them are necessary, the "invariant" condition being a restriction on the sign of the two 
eigenvalues of the Hessian matrix (for example, ()3.2|1 and ()3.4j) imply ()3.3|l ). 

In standard thermodynamics the violation of any of the above inequalities is taken as 
an indication of instability and a phase transition. They are equivalent to the more familiar 
criteria of positivity of the various specific heats and susceptibihties (for example, ()3.2|) 
implies positivity of the specific heat at constant J, etc.) 

3.1 Stability of Black Holes and Non-Extensive Systems 

Note that the above argument relies on the additivity property of the entropy when merging 
subsystems and, more fundamentally, on the very possibility of describing a given system 
as made up of constituent subsystems. Standard thermodynamical stability criteria fail 
in general when applied to black holes because their entropy is non-additive (the BH 
entropy is never an homogeneous first-order function of its variables), and because in GR a 
local definition of mass and other familiar "should-be-extensive" quantities is not possible. 
Therefore, BHs cannot be thought as made up of any constituent subsystems each of them 
endowed with its own thermodynamics, and as such they are non-extensive objects. Let us 
recall that extensivity of the entropy is a postulate in ordinary thermodynamics and 
many general results in thermodynamic theory just fail if this property is not assumed. 

The study of non-extensive thermodynamics has been worked out to some extent, 
though (see e.g. [20] and references therein), and a concrete proposal to avoid the prob- 
lems just mentioned in the BH case has been put forward in a series of articles by several 
authors [121 UHl El ■ These works are based on the so-called "Poincare method of stabil- 
ity" [25 . This method was previously applied to gravitating systems in ^Tj, where it was 
shown to work in the task of showing up the existence of unstable modes of an equilib- 
rium solely from the properties of the equilibrium sequence, without having to solve any 
eigenvalue equation. Formal proofs of the method include those given by Katz |3T] and 
Sorkin [321 E] • A somewhat alternative and more intuitive interpretation has been given 
in [2ni- Let us review these stability analyses before applying them to the BH/BR system. 



As emphasized in one has to consider a single black hole as a whole system. In 
this way, any discussion related to the possibility of dividing it into subsystems or to the 
additivity property of the entropy simply does not take place. Second, one has to put the 
black hole in a given thermodynamic ensemble. The values of the fixed thermodynamic 
quantities that specify the ensemble will be the usual global quantities that one can define 
in GR, thus avoiding any problem related to the non-locality of their definitions. This is 



Setup 
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just what we will be doing if, for example, the chosen ensemble is the micro canonical and 
our starting point is the fundamental relations ()2.6|) and ()2.12j) . 

Considering the BH as a whole system means that, for instance, negative specific heat 
at constant J only indicates an instability of the canonical ensemble: small interchanges of 
energy between the BH and an externa/ thermal bath render the system unstable. However, 
we do not expect such considerations to be related to the issue of classical stability of the 
gravitational solution under small perturbations. 

Equilibrium States, Entropy Maximum Postulate and Fundamental Relations 

The next step is to make a stability analysis based only on maximization of entropy. In 
particular, we do not want to assume additivity of the entropy function. Before proceeding 
to such an analysis, let us first make some definitions and set some notations (we follow 
closely Sorkin [221 ESI)- 

• Let 7W be the space of possible states or configurations of a given system. A given 
configuration is specified by a point X in A^. Let denote the set of independent 
thermodynamic variables that specify a given ensemble, and S be the corresponding 
Massieu function^. Both S and the set of {/i*} are functions on M. which we assume 
to be smooth (for our purposes, having at least continuous second derivatives) except 
maybe in some boundary of A^^. 

• Equilibrium states (stable or not) occur at points in A4 which are extrema of S 
under displacements dX for which d/i* = 0. At any equilibrium state one can define 
conjugate thermodynamic variables such that, for all displacements dX, 

dS = pidfi' , (3.5) 

reflecting the fact that, at equilibrium, the Massieu function only depends on the val- 
ues of the thermodynamic variables {/x*}. Unless the relations dfi^ = are not linearly 
independent, Eq. ()3.5|) defines the quantities {Pi} uniquely. The set of equilibria are 
thus a submanifold A^eq of the configuration space. Curves in this submanifold are 
called "equilibrium sequences", and points in A4eq can be labeled by the correspond- 
ing values of the {/i*}, which are referred to as "control parameters" ''. The points in 
\ A^eq are called off- equilibrium states^. 

'^In the microcanonical ensemble the {/i*} are variables like the energy and the Massieu function is the 
entropy. In other ensembles, the Massieu function is given by the Legendre transform of the entropy. 
^In practice this will mean "extremality" (see below). 

^These names come from catastrophe theory, which is the mathematical framework for the study of the 
phenomena will be interested in. See e.g. 33 

^More properly, one should speak about constrained equilibria, since in principle an entropy function 
is only guaranteed to exist at equilibrium. It is of course difficult to say what "constrained equilibria" 
means in a gravitational setting. In this respect, we will have to assume that thermodynamic functions 
are somehow well defined for off-equilibrium states (see e.g. |34l I35j ). at least in the vicinity of Aloq- 
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• The entropy maximum postulate states that unconstrained locally stable equilibrium 
states take place at points in Aicq where 5 is a local maximum w.r.t. arbitrary 
variations dX in Ai that preserve dfi^ = 0. 

• A fundamental relation for the equilibrium states is an explicit expression for the 



i.e. the integral of ()3.5j) . Given a fundamental relation, the variation of the Massieu 
function as given by Eq. ()3.5p serves to calculate explicit expressions for the conjugate 
thermodynamic variables at equilibrium: 



The above equations constitute the equations of state of the system. 

The entropy maximum postulate is therefore a statement about the behaviour of S along 
off- equilibrium curves in J\4, not about its variations along A^eq- In general, the properties 
of S along A^eq do not tell us anything about the local stability of an equilibrium, and 
expressions like ()3.2|) where S is the equilibrium entropy do not provide information on 
local stability if one cannot invoke additivity of the entropy over constituent subsystems. 
However, the equilibrium entropy as given by ()3.(ij) is often all what we know. In particular, 
it all what we know when considering BHs, because expressions such as fl2.(jj) or ()2.12j) only 
provide values of the entropy function restricted to a very particular region of configuration 
space: namely, the series of (equilibrium) black hole solutions. 

Therefore, questions about local stability based on maximization of entropy (or Massieu 
function) are not questions about the properties of the restricted function Seq- What one 
needs is an extended Massieu function, henceforth denoted by S, that reflects the behaviour 
of the quantity S off the black hole states, at least in their vicinity. However, the form of 
S is not known in general. To have it, one would need a parametric family of spacetimes 
describing deformations of the black hole series, and then to compute all the relevant 
thermodynamic quantities from them. 

The "Turning Point Method" 

The so-called "Poincare" or "turning point" method allows, under certain assumptions, 
to infer some information about the behaviour of S and on stability by using only the 
equilibrium equations of state. In particular, it does not require any knowledge of the 
explicit form of the extended Massieu function. We sketch here the proof originally given by 
Katz in [31] . It has been reformulated for systems like BHs with several control parameters 
in [3ni • Other proofs have been given by Sorkin [321 ^| ■ 

Assume that an expression for the extended Massieu function is given, at least in the 
neighbourhood of Aieq- We can parametrize such a function as: 



function 




(3.6) 




(3.7) 
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where the {X''} denote a set of off-equihbrium variables. Equihbria take place at the values 
X^q = X^q(/x*) which are solutions of: 

a,5.^ = 0. (3,9) 

The equilibrium Massieu function ()3.6p is then given by iSoq(/i*) = 5(Xfq(yU*); /i*). Define 
off-equilibrium conjugate variables (3i, 

^ ^ dS 

A(X'';/x^)^a,5^— . (3.10) 

These give the equilibrium values of Eq. ()3.7p when evaluated at equilibrium, that is, 
I3i{ii^) = Pi{X^^{fi^); fi^). The point represents a stable equilibrium if and only if it is 

a local maximum of S at fixed {/i*}. Therefore, one has to consider the eigenvalues of the 
matrix 

dpJiX^;fi') . (3.11) 

A change of stability along X^q(yU*) takes place when one of these eigenvalues becomes zero 
and changes sign. Assume now that the variables {X''} can be chosen in such a way that 
the above matrix is diagonal. The eigenvalues we are interested in are then: 

A,(/iO = 9^5(X-/iO ^ . (3.12) 

These are called "Poincare coefficients of stability" . An instability appears whenever there 
is some Ap > 0. What we will show now is how it is possible to get information about 
the sign of an eigenvalue near a point where it vanishes without having to compute the 
spectrum of {dpuS)eq. 

Consider the identities provided by equations ()3.9p when the l.h.s. is evaluated at X^^. 
Choose some /i" (for fixed a) and take the partial derivative with respect to it at both sides 
of ()3.9j) . Because of the choice of normal coordinates, one gets the set of relations: 

id,A)eq + XpdaX^^ = 0, (3.13) 

where there is no sum over p in the second term. Now consider Eq. ()3.10p for (with the 
same a as before) also evaluated at X^^. Taking again the derivative w.r.t. /x" gives: 



dp" 



(no sum over a). The important thing is that the l.h.s. contains only the derivative of the 
equilibrium equation of state. Using ()3.13|1 we get: 

|| = 0JSU-5:M, (3,15) 
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Assume for a moment that at a generic equilibrium (c?^iS)eq is finite and that the coefficients 
{dpPa)eq are finite and nonzero (the failure of these assumptions will be discussed below). 
In such a case, the last equation shows that when a change of stability occurs (i.e. when one 
or several eigenvalues Xp become zero and change sign) the plot of the conjugacy diagram 
/3a(/i") along the equilibrium series has a vertical tangent, since dPa/dfi'^ has to diverge 
there — this is a "turning point" (see Fig. |2I). The first consequence of this is that: 

(i) If one single point of an equilibrium sequence is shown to be fully stable, then all 
equilibria in the sequence are fully stable until the first turning point is reached. 

Let us consider the general case in which one can have degeneracy of zeromodes at a given 
equilibrium state, and let p = 1, ■ • • , iV label the eigenvalues becoming zero. Near the 
turning point one has: 



This shows that: 

(ii) At a turning point in the diagram /?a(/i°), the branch with negative slope near the 
turning point is always unstable, 

since a negative slope near the turning point means that at least one term in the r.h.s. 
of ()3.1(ij) has to be negative. 

In a generic case, however, we cannot conclude that the branch with positive slope near 
the turning point is stable: for example, there can always be positive eigenvalues that do 
not change sign and, in such a case, they will never show up in a conjugacy diagram. So 
the positive slope branch is not guaranteed to be fully stable unless we do know that at 
least one point in the branch is in fact fully stable. On the other hand: 

(iii) If the spectrum of zeromodes is nondegenerate, the degree of stability changes by one 
at a turning point. The branch with positive slope near the turning point is always 
more stable than the one with negative slope. 

Here "more stable" means that it has exactly one negative mode less than the negative 
slope branch, since a nondegenerate spectrum means that no more than one eigenvalue is 
changing sign at a time. The positive slope branch will be stable against the particular 
eigenvalue that changes sign at the turning point and, in that sense, it is "more stable" . 

In any case, what one can always see with this method is the existence of instabilities: 
the existence of a turning point from which a negative slope branch emerges necessarily 
implies the existence of at least one unstable mode. Point (ii) above asserts that this conclu- 
sion remains true even without the assumption of a nondegenerate spectrum of zeromodes 
and without any knowledge of stability of any particular point in the sequence IT^ . 




(3.16) 
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Bifurcations 

However, it is still possible to have an eigenvalue Xp that changes sign but whose behaviour 
does not show up as a turning point in the conjugacy diagrams /3a(/i'^). This can only 
happen if, at the point of change of stability, the coefficient {dpj3a)eq also vanishes. Then 
we will not see in general any vertical tangent. It can be shown that this can only happen 
at a bifurcation, i.e. when there is another sequence of equilibria that crosses the considered 
one, the intersection point being the bifurcation. One can indeed prove that P7] : 

(iv) Changes of stability can only occur at turning points or bifurcations. 

The general result is therefore that, in the absence of bifurcations, the stability char- 
acter of a sequence cannot change until the equilibrium curve meets a turning point. In 
particular, stability does not change when the slope of /3a (/^") changes sign when passing 
through zero, thus ruling out the criteria based on the sign of the specific heats (see e.g. 
Eq. ()3.2|) ) for nonextensive systems. This is because, in general, the sign of the slope of 
the diagram Paifi'^) does not say anything about the stability character of any mode. It is 
only in the vicinity of a turning point where both things can be related via Eq. ()3.16|1 . 

Vertical Asymptotes 

One might have wondered why, from the appearance of a zero eigenvalue in ()3.15j) . we 
infer the existence of a vertical tangent, rather than a vertical asymptote. The reason 
is that at a vertical asymptote there are no equilibrium states at all, and therefore they 
signal the (asymptotic) endpoint or "boundary" of the equilibrium sequence. In particular, 
one cannot speak about any "change of stability" there, since if equilibria exist beyond a 
vertical asymptote they belong to a different sequence. 

On the other hand, it seems more reasonable to expect the divergence at an asymptote 
to be a consequence of the non-differentiability of S at the boundary of the sequence, rather 
than to an eigenvalue approaching zero. The 4 dimensional Kerr BH discussed below is 
an example where indeed one can prove that the divergence at a vertical asymptote has to 
be due to some non-differentiability of S. We will consider in detail vertical asymptotes 
in the cases studied below, and we will see that there are important differences w.r.t. the 
behaviour at turning points. 

In the following we will admit non-smoothness of S only at such points. If one does 
not have at least continuous second order derivatives then thermodynamic coefficients do 
not exist, and the Gaussian approximation to the entropy breaks down. 

Examples 

A typical conjugacy diagram is shown in Fig. El The point P is a turning point. The 
upper branch is unstable, while the lower branch can be stable or more stable. The sign 
of the slope also changes at the horizontal tangent at Q, but this has no relation with a 
change of stability. Also, if somehow one knows that there are no bifurcations, then one 
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can conclude that the stabihty of such a curve can change only at P. This means that 
the lower branch would preserve its stability character beyond Q, even if the slope changes 
sign there. 

Let us note that here is a simpler reason to discard the point Q as a point of stability 
change along the equilibrium series. When crossing Q, the control parameter /i" does not 
remain constant, so stability considerations based on maximization of S simply do not 
apply there, since S is maximum at constant /x* (note, on the other hand, that d^"^ = 
at P). This was essentially the reason invoked in pS^ to disregard horizontal tangents in 
conjugacy diagrams. 




p 



Q 



Fi gure 2: Generic plot of a conjugate variable Pa against /i° along an equilibrium sequence. The point 
P is a turning point, where a change of stability can take place. No changes of stability occur at Q, even 
if the slope changes sign there. 

As a concrete example to make the previous arguments explicit, consider the four di- 
mensional Kerr black hole in the microcanonical ensemble . The Massieu function is the 
entropy S and the control parameters are {/x*} = {M, J}. The conjugate thermodynamic 
coordinates are {Pi} = {P,u}, [3 being the inverse temperature and uj = —Q/T. For the 
study of stability we need the plots of /3(M) at constant J and a; (J) at constant M along 
the equilibrium series. Such plots are given in Fig. 01 

The minimum at Md in the plot of (3{M) is exactly the point discussed by Davies [TOj . 
There the specific heat at constant J changes sign through an infinite discontinuity, a 
fact that has been interpreted in the literature as a sign of a phase transition. However, 
this point is a horizontal tangent, and therefore no changes of stability take place at M^. 
Actually, these plots show that there are no turning points at all. Therefore, according 
to this analysis, no changes of stability occur in all the way down from the Schwarzschild 
limit to extremality. There could be bifurcations, but in this case these are ruled out 
due to the uniqueness of time independent BH solutions of the Einstein equations in four 
dimensions. Since we know that one point of the curve (namely, the Schwarzschild limit) 
is stable ^Ei, we can conclude that stability is preserved along the whole curve, at least 
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Figure 3: Conjugacy diagrams for the four dimensional Kerr black hole |15| . 

yS(M) CJ(]) 



M 



Mex, Md 



(a) (3{M) at fixed J. A/cxt stands for the mini- 
mal mass at fixed J (the Kerr bound) . M ^ oo 
is the Schwarzschild limit. Mo is the "Davies' 
point" Jl = (2V3-3)M4 [Xni. 



(b) w(J) at fixed mass. Jext stands for the Kerr 
bound. As in the plot of /3(M), one finds a ver- 
tical asymptote there. 



against the dynamical modes accounted for in this analysis (see the discussion below). The 
non-appearance of instabilities in these diagrams is a result which clearly agrees nicely with 
the fact that Kerr black holes are actually stable jHU]. 

One should consider the origin of the divergence at the vertical asymptote at the Kerr 
bound in view of of Eq. ()3.15|1 . If the divergence at extremality was only due to some Xp 
approaching zero, then that eigenvalue had to be positive (i.e. unstable), since the slope 
is negative. The absence of turning points in these curves would imply that this unstable 
mode never changes sign along the whole series, thus meaning that every point is unstable. 
This is impossible due to the fact that one point of the curve, the Schwarzschild limit, is 
stable. So the divergence at extremality (or, at least, the leading contribution to it) has 
to be due to the non-differentiability of S at the Kerr bound. This makes sense, since 
the extremal black hole has zero temperature and (3 and uj themselves diverge there. So 
even without knowing about the behaviour of 5, one should expect a divergence due to 
the "zero temperature contribution". 

In this case one can say even more about the behaviour of S near the vertical asymptote. 
In principle, both coefficients (9^iS)eq and {dp(3a)eq in fl3.15j) can diverge. However, given 
the fact that we know that all the eigenvalues Xp have to stay negative near extremality, 
the fact that d^/dfi"" is negative there implies that, near extremality: 



{dlS)eq + --- (subleading terms). 



(3.17) 



Note that the above approximation means that the behaviour of the phase diagrams l3a{^Ji°') 
is controlled by completely different physics at a vertical asymptote and at a turning point 
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(compare to ()3.16|) ). This will be important in Section |3] when discussing fluctuations and 
trying to isolate the relevant degrees of freedom both at turning points and at vertical 
asymptotes. 

3.2 Relation to Dynamical Stability 

One should ask about a precise relation between the study of stability based on the Poincare 
method and the dynamical stability under linear perturbations of the black hole solution. 
There are two issues one has to consider at this point: the choice of the ensemble and the 
number and nature of fluctuations. 

Choice of the Ensemble 

We will be interested in the stability of solutions like ()2.1|) and ()2.7|) . These describe isolated 
BHs, for which the mass and angular momentum are conserved quantities that remain fixed 
under Lorentzian time evolution. The appropriate ensemble is thus the microcanonical, 
and the relevant potential is the entropy. Inequivalence of the stability limits of different 
ensembles in self-gravitating systems has been discussed in [TK| ITT)| IT7| EEl HHj and, from 
a more general point of view, in This becomes clear from the precedent discussion 
on the existence of turning points: a different ensemble means a different choice of the 
Massieu function and, typically, a reciprocal choice of the conjugate variables, which may 
in general interchange horizontal tangents ("Davies' points") by vertical ones (turning 
points). This is what would happen if, for example, we decide to put the above discussed 
d = 4 Kerr black hole in the canonical ensemble. In that case stability would be linked to 
the maximization of the "free energy" (p = S — j3M. The conjugacy diagram of interest 
would be M = M(/?), and such a plot has a vertical tangent where (5 = (3{M) happens 
to have a horizontal one. Therefore we see, as mentioned at the beginning of this Section, 
that the point Md only means a change of stability of the black hole when put in contact 
with a thermal bath, but not of an isolated BH. 

Nature of Fluctuations 

As discussed above, the appearance of a single unstable mode at a turning point is enough 
to ensure instability of a black hole solution. However, we are still faced with the ques- 
tion of which dynamical (metric) fluctuations are accounted for in a conjugacy diagram. 
Put differently, we would like to know which deformations of the metric carry physically 
the off-equilibrium fluctuations denoted above by {X''}. In general one will have some 
identification of the kind 

where the {5g\i!v} stand for a suitable set of metric perturbations. This has to be so because 
there are no other possibilities for coordinates in the configuration space^. The problem is 

^For an explicit derivation in specific cases of these kind of identifications from a gravitational path- 
integral point of view see e.g. |40[ I35| and references therein. 
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to identify which particular perturbations are encoded in our analysis. 

This question is important because, in principle, we do not have enough arguments to 
conclude that full Poincare stability (if achieved) is enough to guarantee full dynamical 
stability. First it seems difficult that in a general case the conjugate thermodynamic 
variables could encode all possible physical metric perturbations. In addition, note that 
applicability of the turning point method requires some assumptions on the spectrum of 
(9po-5)eq: namely, a discrete spectrum and the possibility of choosing normal coordinates for 
the off-equilibrium variables. Without any knowledge about the off-equilibrium function S, 
and in the lack of a rigorous proof, one has to admit the possibility that these requirements 
may be, in a generic case, too strong to be satisfied by the whole spectrum of (dpa^S) cq 
under arbitrary linearized perturbations. There could be some perturbations that do not 
satisfy the required conditions for them to be "scanned" with this method, and some of 
those could be unstable modes. 

Perturbations satisfying the above criteria are probably those keeping some degree of 
symmetry^'^. In [T^ IT^ it has been argued that, for rotating BHs, Poincare stability en- 
sures stability against axisymmetric perturbations^^. It would be nice to have an explicit 
proof of this. 

Therefore we see that the turning point method seems reliable when detecting the 
appearance of instabilities. Moreover, one of its main advantages is its extreme simplicity 
and, most importantly, the fact that it can be applied to non-extensive systems such as 
BHs. The drawbacks of this method are, however, the following: 1) Bifurcations do not 
show up in the conjugacy diagrams. 2) Branches with a positive slope at a turning point 
are not ensured to be fully stable unless one can prove stability of at least one point in 
the branch. 3) If it cannot be proven that the spectrum of zeromodes is nondegenerate, 
then one has no control on the relative number of negative modes between two branches 
meeting at a turning point. 4) In principle, there might be dynamical modes that this 
treatment of fluctuations does not encode. 

Before ending this Section, let us mention that this stability analysis only requires 
from the considered system an evolution dictated by the extremization of some potential. 
In our case, such a potential is the BH area, which we know that obeys (classically) an 
extremum law due to Hawking's area theorem. In particular, application of this method 
does not require interpretation of the black hole area or surface gravity as a thermodynamic 
entropy or temperature. 



^"^On physical grounds, it is also unlikely that asymmetric perturbations will extremize the entropy 
(see e.g. @2|)- In such a case, the most natural possibility is that the system meets another branch of 
less-symmetric equilibria (i.e. a bifurcation). 

^^See also |43| for another example concerning the symmetry characterizing the relevant perturbations. 
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3.3 Stability of Black Holes and Black Rings 



Let us now apply these considerations on stability to the five dimensional black hole/black 
ring system. Both for the sake of comparison and for later use, let us turn first to the 
study of the Hessian of the entropy. 



3.3.1 Standard Thermodynamical Stability 

We denote the elements of the Hessian by Hij = d^S/dfi^dfi\ where the variables {/i*} are 
{M, J}. For the rotating BH the fundamental relation that we need is given by ()2.6|) . In 
terms of the variable x defined in ()2.5|1 and the mass parameter /i one finds that (from now 
on we set G = 1): 

23 l-4x' 2V2 1 23v^ 1 

Hmj - 2)3/2 - +^^-2;;7^ + ^(^ 

JJ ^3/2 (1 _ ^2)3/2 ;x3/2 e3/2 2/^3/2 ^1/2^ ^ ^' 

2^ 1 2^ 1 2"^ 1 

where in the r.h.s. we have quoted the divergent part in the expansion around e = 1 — x, 
which we shall need later. The plots of the diagonal elements and the determinant at fixed 
mass are given in Fig. |3] 




Figure 4: Hmm{x), Hjj{x) and AeiH{x) for the rotating BH. The curve changing sign at x = 1/2 is 
that of the element Hmm- Hjj and AeiH are always negative. 



The element Hmm is proportional to minus the inverse specific heat cj at constant J, 
so one can see that cj changes sign from negative to positive values through an infinite 
discontinuity at x = 0.5, exactly as in the four dimensional case. Therefore the same 
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comments about what happens at x = 0.5 apply here (see Section IXTj) . On the other hand 
one would conclude, according to the standard criteria, that the system is unstable for all 
values of the angular momentum, since the Hessian always has one positive eigenvalue. 

As for the BRs, the fundamental relation that we need is given by ()2.12|) . In terms of 
the parameter i^, the elements of the Hessian are given by: 
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(3.19) 
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where one only has to replace z/ by the appropriate value usi^x) or ul{x) in ()2.13|) . Recall 
that z/ = 1/2 gives the ring with lowest angular momentum, z/ = represents the infinite 
angular momentum limit of the large black ring, and z/ = 1 is the x = 1 limit of the small 
black ring. For later use, let us write down the divergent part of the expansions around 
5 = 2v — 1 and e = 1 — z/: 
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Fig. El is a plot of the different quantities at fixed mass for both the small and the 
large black ring. Again, we can see that for both configurations there is always a positive 
eigenvalue of the Hessian: the SBR has one positive eigenvalue, while the LBR has two 
positive eigenvalues for all x. The divergence of the small black ring at x = 1 was expected, 
since this solution approaches the BH solution in the extremal limit (see e.g. Fig.^. Notice 
the infinite discontinuity at x^i^ between the small and large black ring cases, even if both 
solutions coincide and are nonsingular at x^am- 



3.3.2 Turning Points and Stability Analysis 

To perform the Poincare analysis and investigate the existence of turning points we need 
now the plots of f3{M) at fixed angular momentum and uj{J) {u = —Q/T) at fixed mass. 
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Figure 5: Plots of Hmm, Hjj and detH as a function of x for both black rings. 
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X, 



X 



(a) Hmm{x), Hjj{x) and detH{x) for the 
small black ring. All quantities diverge both 
at the lower bound x^i^ and at the upper 
one X = I. 



(b) Hmm{x), Hjj[x) and AetH{x) for the 
large black ring. All quantities diverge at the 
lower bound a;,nin and tend asymptotically to 
zero at a; ^ oo. 



For the black hole and black rings, the conjugate variables are respectively given by: 



All these curves, for the BH and BR phases, are represented in Fig. IHl 

Focusing first on the BH curves (dashed lines) we see that there are no turning points, 
and therefore no changes of stability are shown along the black hole series. Given the fact 
that the five dimensional Schwarzschild solution is stable 03], we can conclude that the 
five dimensional Myers-Perry BH with just one angular momentum is stable, unless: 

a) there are bifurcations, which would imply non-uniqueness of the BH solution other 
than black rings, or 

b) one admits the possibility of the existence of dynamical modes that the turning point 
method does not "see" . 

With respect to point a) above, let us point out that the point x ~ 0.943 in Fig. [T] is not 
a bifurcation, since the BH and the LBR are different spacetimes. One can check that 
the function /?bh = /5bh(^) has a minimum (horizontal tangent) at M = 277rJ^/8 (i.e. 
X = 1/2) which is not shown in the figure^^. This is the "Davies' point" where the specific 

^^The whole BH curves are analogous to those of Fig. |31 




(3.21) 
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Figure 6: The conjugacy diagrams /3(M) and uj{J) for all phases of the BH/BR system along their 
corresponding equilibrium series. The dashed lines represent the BH-curves. We see the existence of a 
turning point in both diagrams at the point where both BR branches meet. The SBR and BH curves have 
a vertical asymptote at extremality. /3bh(-^^) has a minimum (not shown in the plot) at the value of M 
that corresponds to x = 1/2. 
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(a) P{M) at fixed J. Mext stands for the mini- (b) uj{J) at fixed mass, 

mal allowed mass at fixed angular momentum in 
the BH and SBR cases (the "Kerr" or extremal 
bound a; = 1). Mmin is the mass at Xmin. 



heat changes sign but, as aheady discussed, no changes of stabihty occur there. On the 
other hand, the extremal hmit is a vertical asymptote, again like in the four dimensional 
case. The five dimensional case is completely analogous in all respects, and therefore we 
can equally conclude that the divergence at x = 1 comes from the non-differentiability of 
S, and that expressions like ()3.17|1 also apply here. 

Before proceeding to the analysis of the BR curves, let us derive some general properties 
which, assuming no bifurcations, will be of use when analyzing the degeneracy of the 
spectrum of zeromodes. The general expressions of the kind of fl3.15j] are, in this case: 
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These are the elements of the equilibrium Hessian. The off-diagonal term (= Omuj = dj(3) 
can be computed in a similar fashion as we did in Section ()3.1|) to compute the diagonal 
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ones. In the equations above, sum over a denotes sum over never-zero eigenvalues, while 
sum over p denotes the sum over all eigenvalues becoming zero at some point. Note first 
that, near such a point, the diagonal elements of the Hessian have to diverge all with the 
same power: 

Hum ~ Hjj ~ - — h ■ ■ ■ (subleading terms) , (3.23) 

Aj denoting the eigenvalue(s) contributing with the highest power to the sum over p. 
Notice that this need not be the case if the divergence was not due to the appearance of a 
zeromode. On the other hand, the leading contribution to the determinant of the Hessian 
is given by: 

det H ~ h j j _ (Mk(^ j \... (,,,) , (3.24) 

thus meaning that, if there are > 1 zeromodes Ap, then the leading contribution to the 
determinant near a turning point will in general go like: 

detif ~ tV + --- (s-t-), (3-25) 
X1X2 

where Aj represents the eigenvalue with the highest subleading contribution to the diver- 
gence of Ha in ()3.23j) ''^^. If, on the contrary, there is no degeneracy of zeromodes, then one 
will always have the nontrivial property: 

Hmm ~ Hjj ~ Hmj ~ dct ~ ^ + ■ ■ ■ (s.t.) , (3.26) 

since the leading divergence in ()3.24j) will exactly vanish, and the contribution of order 1/Aj^ 
to the off-diagonal term cannot be cancelled. Moreover, in view of ()3.22|) . the numerical 
coefficients h^^^ of the leading contributions to the divergence of the elements Hij have to 
satisfy the exact relation: 

^"-MjJ —"'MM"'.JJ- \0.Z,I j 

This remains true order by order in degree of divergence also in the case of a degenerate 
spectrum in which all zeromodes contribute with a different degree of divergence. 

Let us turn now to the analysis of stability of black rings. From the conjugacy diagrams 
we see that there are two special points: Xmin and x = 1. At Xmin we find a turning point. 
Taking into account the general results of the preceding Section, we therefore conclude 
that the small black ring is always unstable. From the entropy diagram in Fig. ^ we knew 
that the small black ring was globally unstable, since it is the phase with lower entropy. 
Now we see that it is unstable also locally. 
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In the case where Aj is degenerate, the divergence will go like ~ 
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What about the stabihty of the large black ring? First we see, from the behaviour of 
the Hessian near the turning point ()3.20|) . that there is only one divergent piece and, in 
addition, that the necessary conditions ()3.26|) and ()3.27|) for a nondegenerate spectrum 
are satisfied. Therefore this provides strong evidence that the large black ring is more 
stable (in degree one) than the small ring. However, as discussed above, we cannot infer 
from this fact its full dynamical stability. Actually, one expects on physical grounds that 
at least the limit of large x may suffer from a Gregory-Laflamme instability: the large 
angular momentum limit describes a very long and thin ring, and so it may be considered 
as a "black string limit" of the configuration fl^ I45j . If there are GL instabilities in the 
LBR, there are three possibilities to explain the absence of turning points in the conjugacy 
diagrams. 

a) The first one is that these diagrams are "blind" against GL modes. This would be 
the case if, as hinted at in fSlE], the dynamical modes accounted for in this analysis 
are only the ones preserving axisymmetry. Notice that instabilities of the Gregory- 
Laflamme type tend to distort the shape of the string along its axis. In our case, this 
means deformations along the within the plane of rotation, and therefore they are 
not axisymmetric perturbations. 

b) The second possibility is that the GL instability is present for all values of x. In this 
case the LBR would be always unstable but there would be no turning points. 

c) The last possibility is that the emergence of GL modes actually appears as a bi- 
furcation. This would imply the existence of a new branch of solutions beyond 
some angular momentum per unit mass, most probably describing non-uniform black 
rings |24| EH] . We consider this latter possibility an interesting one which deserves 
further investigation. 

A further heuristic argument, due to Emparan and Reall, in favour of another instability 
of the large and small black rings is the following. As discussed in jlH], one expects the 
limit X = Xmin to be unstable: one simply has to throw matter with no angular momentum 
into the ring with minimal x. This would decrease the value of the angular momentum per 
unit mass of the final state and, given the absence of black ring states at x ^ -fmin? such 
a process would destabilize the system. This could happen if possibilities analogous to a) 
or b) above apply to such an instability. 

Let us comment now on the origin of the divergence at the vertical asymptote that the 
SBR curves exhibit at x = 1. As discussed for the black hole case, it seems most likely that 
such a divergence appears due to the fact that P and lo themselves diverge there because of 
the zero temperature at x = 1. In the SBR case, however, one has to admit the possibility 
that the divergence can be due to an unstable mode that approaches zero at extremality 
since, contrary to the BH, the SBR has negative modes for all x. We want to argue here 
that this does not seem the case. The BH and SBR solutions approach each other as they 
get close to X = 1, so one should expect the physical origin of this divergence to be the same 
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in both cases. In fact, by looking at the leading contributions (~ e^^/^) to the divergence 
in the expansion of the SBR Hessian ()3.20|) at x = 1, we see that these exactly coincide 
with those of the BH ()3.18|1 up to an overall constant factor. So we find evidence that the 
divergence at a; = 1, or at least the leading contribution ~ e~'^/^, is due to the fact that 
the temperature vanishes there. Therefore, as in the four dimensional case, we will write: 

Hij ^ {dijS)eq + ■ ■ ■ (subleading terms) , (3.28) 

both for the BH and the SBR near extremality (we are assuming that this expression also 
holds for the off diagonal term of ()3.22|) ). This will be of importance in the next Section, 
when we will consider the behaviour of fluctuations. 

Summarizing, our analysis proves local instability of the SBR and presents strong ev- 
idence that the LBR is a "more stable" phase, even though it does not guarantee its full 
dynamical stability. In fact, one expects the LBR to suffer from instabilities which do not 
show up in the conjugacy diagrams, the possible reasons for this being the ones discussed 
above. On the other hand, these conjugacy diagrams do not show any special property nei- 
ther at a; = 1/2 ("Davies' point") nor at x = 2-\/2/3 (where the LBR becomes entropically 
favoured, see Fig. [TJ. 

4 Fluctuations, Critical Behaviour and 
Thermodynamic Geometry 

When an unstable point is reached one expects fluctuations to diverge, as it also happens for 
instance at the critical point of a second order phase transition. However, both phenomena 
are not the same. Let us next elaborate on this issue in some detail. 

Before proceeding we want to make clear what we mean by "fluctuations" in the con- 
text of the off-equilibrium formalism introduced above. What we want to consider here are 
general motions in configuration space, i.e. motions parametrized by {AX'', A/i*}. Note 
that for the study of stability we needed to work at fixed {/x*}, as required by entropy max- 
imization. Here, however, we will not be interested in stability considerations, but rather 
in the changes induced in the system as we move in all possible directions of configuration 
space. This is why now we will consider variations like {A/z*} as well. 

The physical significance of the variations in the {X''}- and in the {/i*}-directions is, 
however, different. The off-equilibrium variables {X''} should be considered as "dynami- 
cal", in the sense that they obey a first order equilibrium equation {dpS = 0). Therefore, 
the variations {AX''} can be thought as spontaneous fluctuations of the system. On the 
other hand, the set of variables {/i*} are our control parameters (in the microcanonical 
ensemble, the mass and the angular momentum) and, as such, they do not fluctuate spon- 
taneously. In any case, we can always consider the changes induced in the system as we 
change these control parameters ( "adiabatically" , say), in an attempt to map out the be- 
haviour of the system (or better, the thermodynamic potential) in configuration space at 
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the vicinity of equilibrium^^. At this point, let us recall that all our considerations are 
purely classical, and therefore we are ignoring any quantum effects like Hawking radiation 
(in such a case, spontaneous variations in mass and angular momentum would be of course 
physical). 



4.1 Fluctuations at Turning Points and Vertical Asymptotes 

Let us fisrt consider fluctuations around equilibria near a turning point, following in part 
the discussions in jlOl EEl ES] • If the eigenvalue becoming zero there is nondegenerate, 
then the fluctuations will be basically driven by fluctuations in the corresponding mode 
X^, since this is the direction along which the Massieu potential becomes fiat, the restoring 
forces are small and then AX^ can be large. This will induce off-equilibrium fluctuations 
in the conjugate variables which will be well approximated by: 

SPi ~ {d-A)e^AX'' (4.1) 

(unless (9j/3j)eq is vanishingly small, which can only happen at a bifurcation — see Sec- 
tioning. This approximation remains valid also if there are other zeromodes, but still the 
one approaching faster to zero is nondegenerate. 

The probability for the system to be at X = Xeq -|- AX^ is proportional to the expo- 
nential of the Massieu function evaluated at X |^ 14 7j . which will be well approximated 
by its second order variation along that direction: 

P(Xeq + AX~') ~ exp + ^XiiAX'^j . (4.2) 

The last two expressions together with Eq. ()3.16|) show that, near a turning point which is 
not a bifurcation and where Aj is nondegenerate, the probability density for a fluctuation 
6Pi is given by 

1 mr 



P{5l3i) ~ exp 



2 dl3i/dii' 



(4.3) 



Calculating the normalization factor and evaluating the mean square deviation (from the 
side with Aj < 0) is straightforward, and gives the result [101 EEl ES] : 

((^A)^) - . (4.4) 

The same result was obtained in [SEl ^| via somewhat different considerations. There the 
were taken from the start as good coordinates to parametrize off-equilibrium states, 
thus casting all the discussion on off-equilibrium fluctuations directly in terms of them 
(without passing through the {X^} coordinates used here). In such a case, the above 

^^Note that variations A^' will not in general run along the equilibrium series unless they are correlated 
with the AXP as dictated by the equilibrium equation X^^ — Xp{ii^). 
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formula, which here is obtained only as an approximation near a turning point, turns out 
to be valid everywhere^^. However, in general one will need a different and more complete 
parametrization to account for all off-equilibrium states. In such a case the above equality 
is ensured to hold only near a turning point. 

So we see that, when the system approaches a turning point from the stable side, 
the quadratic fluctuations of the conjugate variables diverge. In the BH/BR system this 
happens at Xmin, where we also expect the condition of Aj being nondegenerate to be 
satisfied (see last Section). In such a case we can write, near Xmin (see ()3.22|) ): 

-"MM ~ ^ , -n JJ ~ ^ , nMJ ~ ^ ■ (,4.5j 

A| Ai Aj 

Let US recall that the dominant term in the expansion of the BR Hessian fj3.20|) exactly 
obeys ()3.27|1 . i.e. the numerical coefficients are in perfect agreement with the approxima- 
tions above. Calculating the ratio Sp/Su from ()4.H) allows to write, from ()4.4|) and Hmj 
above, the off-diagonal correlations (5/5 6u) . We thus obtain the general expression: 

{5Pi6^j) ^+H,j (4.6) 

for correlations near a turning point. The expression above is then telling us how fluctua- 
tions of the LBR diverge as the system approaches Xmin- 

Let us turn now to the fluctuations in the BH and the SBR cases near the vertical 
asymptote at a; = 1. Here the phenomenon is very different, since to start with is not 
related to a change of stability. We saw in the last Section that the leading divergence of 
the elements of the Hessian can be explained by the divergence of the elements {dijS)eq- 
Eq. ()3.28j) implies that these coefficients are the ones that dominate also the series expan- 
sion of the fluctuations 6(3 and 6uj (since it implies that {dijS)eq ^ {dpPi)eq {dpPj)eq) and, 
in addition, that they are well approximated by the equilibrium Hessian. Therefore we can 
write: 

^ Hi.Afi^ . (4.7) 

So we see that, contrary to the previous case (see Eq. (HHl)), the variations 6(3i are now 
dominated by the variations in the variables {/i*}, which are small but still induce the main 
contribution to because of the large coefficient. Variations in the {/i*}-directions will 
be governed by a probability density which, using ()3.28|1 . can be written as approximately 
proportional to 



P{fx' + Afx') ~ exp 



^^This can be seen from the general formalism by substituting in l|4.2|l (which in this case holds for all 
modes, which are precisely the S/Si) the exact value of {d?S)cq as obtained from Eq. (|3.15() and taking each 
Pi as independent from the {/i*} and the rest of the f3j, j ^ i). 
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Expressions ()4.7|) and ()4.8|) formally coincide with those used in the usual treatments of 
fluctuations in thermodynamic fluctuation theory^^ jlTj. The computation of correlations 
gives the result: 

(<5A 5%) ^ -Hij + PiPj ^ -Hij (4.9) 

(note the minus in the final expression sign w.r.t. ()4.fij) ). The computation of these fluc- 
tuations requires some care, however, due to the presence of a positive eigenvalue in the 
Hessian matrix. Its derivation is outlined in the Appendix. The first equality is the exact 
result as computed from the (approximate) probability distribution ()4.8|) . The additional 
contribution to the Hessian is due to the presence of the linear term in the probability 
distribution. However, such a contribution can be neglected near extremality, since it is 
subleading as e ^ (see Eqns. flTT^ . (IT^ and (IT^ ). 

Summarizing, we see that both near extremality and near the turning point fluctua- 
tions diverge, and that at those points correlations can be computed from the equilibrium 
Hessian. 

Another situation where fluctuations typically diverge is at a critical point. However, as 
pointed out in ^HlEl; this is in general a different phenomenon from a change of stability. 
A critical point is the endpoint of a "coexistence curve" , characterized by the existence of 
two competing local maxima of the corresponding Massieu function, describing for example 
the liquid and gas phases of a fluid. At the critical point both maxima coalesce, which 
entails a locally flat potential and an absence of restoring forces. This is what makes 
fluctuations become large near the critical point. On the other hand, the divergence of 
fluctuations related to a change of stability is in general due to a different mechanism: the 
flattening of the potential because the system enters into a series of local minima (unstable 
equilibria) from a series of local maxima (stable equilibria). 

Next we would like to investigate these phenomena in the BH/BR system. To this 
aim, further information on the critical behaviour of a system can be obtained from the 
geometry of the thermodynamic state space, which is what we shall introduce next. 

Before ending this Section the following comment is in order. Note that, at least in the 
SBR case (and probably in the LBR also), we are considering fluctuations of an unstable 
system (i.e. with some Ap > 0) . Therefore, to make sense of all approximations above, we 
have been forced to suppress by hand all variations along the unstable {X^}-directions, 
which will be of course the dynamically relevant ones. All our considerations are thus 
purely formal, in the sense that we are restricting ourselves to study these configurations 
along the directions in configuration space in which no unstable modes are present. It is in 
this sense in which this analysis on fluctuations (and the subsequent discussion on a phase 
transition) has to be understood. 

^®Up to the linear term in the probabihty distribution, whose effect just amounts to a shift of the mean 
values (<5/3i) (see Appendix). As shown below, the contribution of this term to the mean square deviations is 
subleading. The linear term is nonvanishing here because the {/i*} do not obey any "equilibrium equation" 
(i.e. the vanishing of a first order variation) — only the {X^} variables do. 
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4.2 Critical Phenomena and the Geometry of the 
Thermodynamic State Space 

It has been observed by several authors in different contexts that the entropy function 
induces a natural metric on the thermodynamic state space, and that the geometric in- 
variants constructed out of it may provide interesting information about the phases of the 
model under consideration. This formalism was pioneered by Ruppeiner, and we refer 
the reader to |2Zj for a very detailed review and various examples. Ruppeiner formalism 
relies on the usual thermodynamic properties of extensive systems. Here we want to argue 
that, modulo some careful reinterpretations, it can be applied also to BHs near turning 
points and vertical asymptotes. Previous works where Ruppeiner and related thermody- 
namic geometries have been applied to black hole physics include jlHj, but the issue of the 
non-extensivity of BHs has not been discussed there. 

Let us first review how this formalism can be motivated in ordinary (extensive) ther- 
modynamics j2ZI- Let us consider the microcanonical ensemble, in which one starts with 
a homogeneous, isolated system in the thermodynamic limit. The purpose is to investi- 
gate thermodynamic fluctuations in a given open subsystem (to be thought as a "small" 
subsystem) of fixed volume V. The microcanonical variables {/i*} are fixed within the 
total system, but the ones specifying the state of the subsystem {/iy} can fluctuate around 
their equilibrium values {/^Veq}; which will be set by those of the environment. All off- 
equilibrium states of the subsystem are assumed to be uniquely specified by the values of 
the {nir}. 

According to the basic assumption of thermodynamic fluctuation theory, the probability 
density Py of finding the subsystem in a state characterized by values {/iy} of the remaining 
variables of state (all except the volume, which is fixed) is given by: 



where the S'tot is the entropy of the total system and /i^q are the average (equilibrium) 
values of the environment. The function S'tot (/^V ' /^eq) has, as a function of {/iy}, a local 
maximum when the densities of the subsystem equal those of the environment, i.e. at 
fiy = /^Veq ~ ^/^eq/^ot- The abovc formula describes fluctuations of the subsystem around 
this equilibrium value. If one now assumes extensivity, i.e. the fact that the entropy of the 
total system Stot is the sum of that of the subsystem Sy plus that of the environment, and 
expands around the equilibrium point up to second order, one gets jl?! : 



where (n — 1) is the number of independent fluctuating variables once the volume has been 
fixed. In the expression above gij is defined as: 




(4.10) 




(4.11) 



(4.12) 



'V 
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and A/iy = /^y — /^Veq- S^'^ '^^^^ repeated indices (i, j = 1, 2, ■ ■ ■ , n — 1) is understood in 
the argument of the exponential. In standard treatments the formula above is written in 
terms of the entropy and other quantities per volume (i.e. Sy = Vs, etc). Here we prefer 
to leave it in terms of the absolute quantities of the subsystem itself to make closer contact 
with the black hole case. 

So far this is standard thermodynamic fluctuation theory [37]. Ruppeiner proposed to 
interpret the Hessian of the entropy metric in the thermodynamic state space. This 

metric is normally known in the statistical mechanics literature as the "Ruppeiner metric" , 
and it is supposed to govern the fluctuations of the system. The meaning of distance as 
defined by this metric is the following: the less the probability of a fluctuation between 
two states, the further apart they are w.r.t. the geodesic distance defined by this metric. 
In other words, the metric is an indicator of the "closeness" of probability distributions. 
A straightforward computation based on the Gaussian distribution shows that the average 
values of fluctuations vanish: 

(A^V) = 0, (4.13) 
and that the main square deviations are given by the contravariant metric tensor: 

(A/iVA/i{.) = (7^^(/ieq) , (4.14) 

which becomes an alternative definition of the thermodynamic metric. 

Ruppeiner theory is covariant, in the sense that the line element built from the metric 
defined in fl4.12j) is taken to be an invariant quantity. This prescription serves to tell us 
what to do if the chosen coordinates to describe fluctuations of the subsystem are not the 
usual extensive quantities. If for example we choose to work with the conjugate variables 
in the entropy representation 

^v^ = ^, (4.15) 

the line element is written as: 

ds^ = gijifiv) dnydnl. = -d(3vidjjiy = J d(3vidf3vj , (4.16) 

OlJViOpvj 

where the new Massieu function cj) is the Legendre transform of the entropy w.r.t. the 
variables {(3vi}'- 

cPiPv) = Sv-l^'vPvi. (4.17) 
Equation ()4.16|) gives the expression for the metric in the new coordinates^^. 

Expressions like ()4.13j] and ()4.14j] regarding fluctuations hold in any chosen coordination 
of the state space. 

^^Our use of upper and lower indices could be confusing here. Notice however that a change of variables to 
the makes them become the new genuine contravariant coordinates. Expression (|4.17|) gives therefore 
the covariant components of the metric in the new coordinates. 
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Thermodynamic Curvature 

A very special quantity in Ruppeiner theory is the curvature computed from the thermo- 
dynamic metric. For a two dimensional system (n — 1 = 2), as will be later our case of 
interest, the curvature scalar contains all nontrivial information about the metric itself (in 
higher dimensional cases one may need other components of the curvature tensor). 

We will not review here the derivations leading to the meaning of the thermodynamic 
curvature scalar. We will simply recall its main properties, referring the interested reader 
to the Ref. j2Zj. First, computation of the scalar curvature for known statistical models 
provides indication that it is a measure of interactions. For example, it is zero for the ideal 
gas, but the thermodynamic metric turns out to be curved in the Van der Waals case. But, 
most importantly for our purposes, the curvature diverges at a critical point, and it does 
so in the same way as the correlation volume i.e. 

i?crit~e', (4.19) 

where ^ is the correlation length and d the effective spatial dimension of the system. The 
reason for this to be so can be put a bit more precisely: calculations show |2ZI that the scalar 
curvature measures the smallest volume where one can describe the considered subsystem 
as surrounded by a uniform environment. At criticality, one expects such a volume to scale 
as 

Applicability to Black Holes 

So far we have sketched the main features of Ruppeiner formalism as they can be motivated 
from ordinary, extensive thermodynamics. We deduced the notion of a metric in the 
thermodynamic state space by considering a subsystem of a much bigger extensive system 
in the microcanonical ensemble. Such a total system effectively acts as a reservoir for the 
subsystem, of fixed volume V. A careful derivation of the theory just presented requires 
extensivity of the entropy function. We have in fact invoked extensivity to arrive to the 
first place where the thermodynamic metric appears, Eq. 1)4.1111 . How can all this be 
conceivably applied to non-extensive systems such as black holes? 

We will try to argue now that this is in fact the case. All considerations exposed so far 
come from the Gaussian expansion 1)4.111) . In statistical mechanics, extensivity is required 
to infer the fundamental relation Sy from the knowledge of S'tot, which is what is usually 
known. It is also needed to "split" the degrees of freedom of the subsystem from those of 
the reservoir. But note that, in practice, this splitting only serves to fix the equilibrium 
values of the subsystem. On the other hand, since we are fixing one scale, the function 
Sy will not be an homogeneous function of its variables. In practice, what this means 
is that we will be working with an entropy function Sy which is not extensive within the 
subsystem. The latter and its fluctuations are the object of physical interest, and the 
formalism described above does not meet in principle any difficulty as long as: 1) the (non- 
homogeneous) function Sy is known, independently of its relation to S'tot, and 2) we have an 
expression for the relevant off-equilibrium variations of the subsystem around its (known) 
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equilibrium state. This is precisely the situation we found for black holes near turning 
points and vertical asymptotes, which therefore have to be identified with "Ruppeiner's 
subsystem" . 

In Ruppeiner theory, general off-equilibrium states of the subsystem are described ex- 
actly in Eq. ()4.11|1 by the values of the {/iy}. That is, all possible off-equilibrium states of 
the subsystem are supposed to be well described by the thermodynamic limit. In the the 
BH case, however, we do not have in general any control on off-equilibrium fluctuations, 
since they are given by the completely unknown function S{Xf\^^). The analog of the 
Ruppeiner metric (i.e. the coefficients of the quadratic expansion of the entropy) will be 
given, in a general case, by the many-dimensional metric with elements: 

9ij = -{^ijS)oc^ , 9ip = -{^ipS)oc^ , gpp = -Xp ■ (4.20) 

However, near a turning point or a vertical asymptote, the expressions ()4.2j) and ()4.8|) 
provide a good description of off-equilibrium fluctuations. Moreover, the quadratic ap- 
proximation to the entropy is computable, and the "effective" elements of the Ruppeiner 
metric are just a few and can be computed from the equilibrium Hessian. In fact, one 
only has to use the prescription ()4.14j) for the inverse metric as read from equations ()4.6|) 
and (Oil . 



4.3 Thermodynamic Curvature for the BH and BR Spacetimes 



As explained above, the elements of the Ruppeiner metric can be computed from the 
expressions for quadratic fluctuations. Near x^i^ and x = 1 we will write, according 
to (gH) and 

g'^ = = ±H,^ , (4.21) 

where the plus sign stands for the LBR at x^in and the minus sign for the SBR and BH 
at X = 1. Let us compute now the curvature scalar from these metric elements. For 
convenience, we give below the expressions as computed from g^^ = +Hij, and we will 
provide their values for all x. The physically meaningful limits are however those at Xmin 
and X = 1 with the appropriate sign. The results for the rotating black hole and black 
rings are, respectively: 



R 



1 



BH 



R 



BR 



2v^ v^(z/2 + 2z/-2) 



(4.22) 



^2^3/2 



{2u-iyvT 



V 



These quantities are plotted in Fig. We see that divergences of the curvature (i.e. 
indicators of critical behaviour according to Ruppeiner theory) are found at x^m and 
x = 1, in agreement with our previous analysis on stability and fluctuations. 

""^^For simplicity, we are neglecting the subleading contribution /3i/3j in (|4.9|l . It can be checked that this 
does not change the power behaviour of the divergence of the thermodynamic curvature at extremality, 
which is what we will be interested in. 
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Figure 7: Thermodynamic curvature scalar for the BH (dashed hue), SBR and LBR phases. The 
physicaUy relevant regions are those near Xmin and x = 1. The curvature is divergent at those points. 



5 Scaling Laws in the BH/BR System 

Ruppeiner analysis confirms that the points Xmin and x = 1 are candidates of critical 
points. However, divergences in the thermodynamic curvature have been found to happen 
also when the system enters into an unphysical region. In j^H], for instance, computation 
of the thermodynamic scalar for the Van der Waals model showed that it diverges also at 
the so-called spinoidal curve. In our case, both at Xmin and at x = 1 the BH/BR system 
enters in an unphysical region (namely, naked singularities). One can also argue that the 
origin of these divergences could just be the divergent nature of fluctuations at those points 
(which, as discussed in Section [4.11 do not necessarily mean critical behaviour). In order 
to clarify this issue, we examine next the scaling relations between the different critical 
exponents both at Xmin and x = 1. 



5.1 Critical Exponents and Order Parameter 

We shall first define the appropriate susceptibilities and critical exponents that are suited 
to the microcanonical ensemble. Following [18j we define them as follows: 



—a — 



(5.1) 



These are the natural generalizations of the expressions familiar from other ensembles (for 
example, in the grand canonical ensemble the thermodynamic potential is the Gibbs free 
energy, and the specific heat is given by its second derivative with respect to the tempera- 
ture, etc). Also, in the expressions above the e parameters are the suitable generalizations 
of what e.g. would be the "reduced temperature" (T — Tcrit)/TCTit in the usual treatments 
(i.e. (grand) canonical ensembles) of more familiar examples like liquid/gas transitions, etc. 
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We want to study the scaling relations at the two points of the system where fluctuations 
diverge, i.e. at Xmin and at extremality, x = 1. At each of these, we define the parameters 
eM and ej to be: 



At 



Cm 



M-M. 



At X = 1: 



Jm\r 



ext 



(5.2) 



where 



J rr 



7rj2 
G 



2 \ 1/3 



ext 



/277rJ 
V 32G 



J ext 
2\ 1/3 



CMS 



TT 



32GM3 
277r 



(5.3) 



Other critical exponents tell us about the behavior of the order parameter of the tran- 
sition considered. In general, there is no obvious choice for an order parameter. In ^H] it 
was proposed to choose, as an order parameter for BH phase transitions, the difference of 
the "intensive" (conjugate) variables between the two phases. In our case this means: 




(5.4) 



These parameters go to zero at the (presumed) transition point, as they should. For these 
order parameters, the natural generalizations of the critical exponents (3 and 6 familiar 
from statistical mechanics is given by: 



V 



-M 



'J 



(5.5) 



With these definitions, from equations ()3.18|1 and ()3.19|1 the critical exponents are readily 
computed to be: 



At X = 1: 



Small Black Ring : a 

Large Black Ring : a 

(3=1/6= 1/2. 
Black Hole : 

Small Black Ring 

(3 = 1/6 = -1/2. 



^ = 1 
= 7 



a = 1/2, 
1/2, 



Q/ = (^ = 7 = o- = 3/2, 

a; = y9 = 7 = cr = 3/2. 



(5.6) 



Incidentally, let us notice that these critical exponents at extremality precisely coincide 
with those found in ^Tj ^j, where four dimensional black holes were considered. In 
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that case, the two proposed phases were the inner and outer horizons, the black hole 
solution being thought then as a "coexistence curve" . In the five dimensional case we are 
considering, our candidates for the two different phases are instead the black ring and the 
black hole. 




5.2 Scaling Laws and Correlation Length in the BH/BR System 

With these critical exponents at our disposal, let us check now if the scahng laws are obeyed. 
First of all, note that a nontrivial result is the equality of the critical exponents a, y?, 7 and 
a at both phases of each transition point considered (something which, in principle, need 
not be so). As for the scaling relations involving them and the other critical exponents 
defined so far, we quote here their expressions. They are given by [TTj : 

= 2, 

= 7, (5.7) 
= a . 

We can verify that these scalings are obeyed by the BH/BR system, both at x = Xmin and 
qXi X —— X . 

However, according to the general theory of critical phenomena, further scaling laws 
involving additional critical exponents have to be satisfied near a critical point |17]. The 
remaining critical exponents are those related to the behaviour of the two-point correlation 
function and the correlation length. The two-point correlation function G{r) introduces a 
new critical exponent ( defined by: 

exp(-r/0 

{d is again the spatial dimension of the system). On the other hand, the divergence of the 
correlation length ^ at the critical point is expressed in terms of the critical exponents 11 
and u as: 

e ~ ~ e-/ . (5.9) 
At a critical point these exponents obey the scaling laws given by |47j : 



(5.10) 



To verify these relations we need to compute the correlation length and the two-point 
correlator. However, there is no obvious way (not even an obvious definition) to compute 
these quantities for a black hole geometry. In any case, Eq. ()4.19p provides us with an 



u{2 - 


= 7, 


2-a 


= ud , 




= u . 
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explicit algorithm to compute the quantities i/d and jj^d if applicability of Ruppeiner theory- 
is assumed. One obtains: 



Relations ()5.1()j) are obeyed at extremality for any effective dimension d, but are not sa- 
tisfied at Xmin- 

We therefore see that x = 1 satisfies additonal scaling laws not obeyed at x^i^. In 
this respect, we can conclude that the behaviour of the BH/SBR system near extremality 
enjoys properties which, formally, are analogous to those of a critical point of a second 
order phase transition in which the SBR is always the "disfavoured" or metastable phase. 
As explained at the end of Section 14. ![ this statement has to be understood as describing 
the properties of the system only along the stable directions in configuration space. 

On the other hand, at Xmin "we do not find an analogous behaviour of the large and 
small ring phases. In spite of the divergent behaviour of fluctuations, we see that the 
scalings related to the critical exponents of the correlation length are not satisfied there. 
This is in agreement with our expectations about the origin of divergent fluctuations at 
Xmin where, as we saw, what we have is a change of stability. As discussed in Section Wl\ 
this phenomenon is very different from that of a critical point. 

6 Conclusions and Outlook 

We have addressed the issues of stability, divergence of fluctuations and critical phenomena 
in the BH/BR system. Our results give a consistent and unified picture of the properties 
of the phase diagram of black holes and black rings. The two main tools that we have used 
are the Poincare method of stability and the geometry of the thermodynamic state space. 
From a quite general point of view, we have put special emphasis on the consequences of 
the non-extensivity of BH thermodynamics, and we argued on the suitability of the meth- 
ods used here for the study of BHs. The main difference between the approach followed in 
this paper and others in the BH literature is that ours is based crucially on the study of the 
off-equilibrium function S'(X^;/i*), and not on the properties of the equilibrium entropy 
S'eq(/i*). The main advantage of the Poincare method is precisely that, just by plotting 
the right phase diagrams of an equilibrium sequence, it can be used to infer some relevant 
properties about S{XP] /x*), even if such a function is completely unknown. 

Let us he summarize the main results concerning the BH/BR system. First, the small 
black ring has been shown to be locally unstable. Concerning the large black ring, we 
found strong evidence that this is a more stable configuration, i.e. one with less unstable 
modes. It is nevertheless true that this does not imply its full dynamical stability and, in 
fact, one expects such a configuration to be unstable. The only place in the phase diagram 
of Fig n where we see a change of stability is at Xmin, where both BR branches meet. 





ud = iJ,d = 1 , 



(5.11) 
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We are able to evaluate thermodynamic fluctuations and compute the second moments 
of correlations near Xmin and x = 1. In this respect, we have discussed in detail the 
differences between turning points and vertical asymptotes. We saw in Section f4.1l that, 
if one considers fluctuations only along the stable directions, near a turning point the 
fluctuations of the conjugate variables are dominated by the less stable mode X^, i.e. 

due to the fact that AX^ is large. On the contrary, near a vertical asymptote the fluctua- 
tions 6f3i are dominated by those in the usual thermodynamic variables, i.e. 

due to the fact that {djj3i) is divergent. This is likely to be a more general phenomenon. 

Finally, the computation of the various critical exponents allows to formally interpret 
X = 1 as a critical point along the stable directions in configuration space, while at Xmin 
the divergence of fluctuations seems to be related only to a change of stability. This fits 
with our expectations from the stability diagrams in Fig. IHl since at a turning point we 
only expect a change in stability. Evidence of a critical point at a vertical asymptote has 
also been found in [13 CH] • 

The extremal configurations at x = 1 are naked singularities. However, the general 
properties of a system at criticality concern its behaviour near the critical point, and 
our configurations are nonsingular there. A critical point indicates the endpoint of a 
"coexistence curve" , along which a first order phase transition between the two competing 
phases can take place. It is true, however, that in the BH/SBR system such a first order 
phase transition would never be physically relevant in the dynamics of these spacetimes: 
not only because the SBR is locally unstable, but also because this phase always has lower 
entropy than the BH. Also, we do not know about any phases beyond the "critical point" . 
In any case, the fact that the scaling relations are satisfied looks nontrivial. It would be 
interesting to understand better the physical relevance of this property. 

On the other hand, the absence of any other points of change of stability (in addition to 
the general arguments to rule out stability considerations based on the sign of the specific 
heat) makes it very unlikely to interpret the change of sign in the specific heat in the BH 
branch as a phase transition. Moreover, at x = 1/2 it is not clear at all between which two 
phases there could be a transition. 

Finally let us comment on the other special point in the phase diagram of Fig.^ namely, 
X = 2v^/3. Our results do not indicate so far the existence of any special properties there, 
neither in the BH nor in the LBR branch. Given this fact, we therefore conclude that if 
both branches are stable around this point, the only possibility we are left with is that 
of a first order phase transition between both. We find no indications at all to interpret 
X = 2^/2/3 as a critical point of a second order phase transition between the BH and the 
LBR. 
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From the point of view of the apphcabihty of the techniques used in this paper, it 
would be nice to understand in a more precise way the relation between the off-equilibrium 
fluctuations {X''} considered here and the dynamical modes {Sg^^} of metric perturbations. 
Given the simplicity of the turning point method, it would be very interesting to find an 
explicit and precise connection, and see which dynamical modes these conjugacy diagrams 
are actually probing. Even if, in general, one cannot expect to recover full information on 
dynamical stability from the Poincare method, it seems clear that this method looks reliable 
at least to detect instabilities. By this we mean that Poincare instability (a negative slope 
near a turning point) necessarily implies the existence of at least one unstable dynamical 
mode (while the opposite is not necessarily true). 

In this respect, a natural question to ask is about the relation (if any) of the Poincare 
method used here and the Gubser-Mitra conjecture ^Hl- The latter relates, for black 
brane solutions, standard local thermodynamical instability to the presence of dynamical 
instabilities in the same sense as discussed here^^. We have argued against the applicability 
to BHs of the criterion of stability based on the sign of the various specific heats. Note, 
however, that the GM conjecture does not include compact horizons, just black branes 
with some noncompact translational symmetry. This means that, along the noncompact 
directions, one has extensivity of the mass per unit length and additivity of the entropy 
per unit length. Therefore standard thermodynamic criteria should apply in these cases, 
at least for perturbations generating lumps along the uniform directions. Note that these 
kind of GL instabilities are in perfect qualitative agreement with the simple example that 
we discussed in the introduction of Section El (a negative Hessian means that system is 
unstable against redistribution of mass and tends to become non- uniform) . It would be 
extremely interesting to explore this issue in more detail and see, in particular, if the 
appearance of this kind of instabilities is related to bifurcations of equilibria. 

Concerning the use of the thermodynamic metric and its relation to BH physics, it 
would be desirable to have a precise definition and understanding of the correlation length 
for a BH geometry. What we have used here is the divergence of the thermodynamic 
curvature. This is probably measuring the typical size of fluctuations in the BH geometry 
as compared to the size of the BH itself. It would be nice to work this out in an explicit way. 

The techniques we have been using are completely general, and obvious extensions 
of this work include their application to other gravitational phase transitions of interest. 
Work on the directions proposed here is in progress. 
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A Fluctuations Near Extremality 

Let us sketch here the derivation of equation ()4.9p for fluctuations near extremaility. In 
ordinary systems with positive specific heat, the probability distribution just reduces to an 
ordinary Gaussian integration due to the fact that the eigenvalues of the Hessian matrix 
are all negative. Here, however, we always have one positive eigenvalue, which makes the 
normalization constant of the probability distribution strictly infinite if one is to integrate 
fluctuations over the whole real line. Therefore, in evaluating the expressions ()4.9|) . we 
need to impose a cutoff in phase space, i.e. 



(5A5/3,-)(Am,Aj) = ^^^^^-^ I d{AM) I d{AJ) 6l3,6(3j exp{S), 

I-Am J-a 



N{Am,Aj) 



pAM pAj (^•^) 

N{Km.Aj) = / d{AM) d(AJ)exp(5), 

J -Am J-Aj 

where S stands for the quadratic approximation to the entropy near the vertical asymptote, 
i.e. Eq. ()4.8p . and the variations SPi are given by Eq. ()4.7p . 

Before proceeding, let us note that a cutoff is always necessary: if one is to make 
sense of any truncated series expansion in A/i* of the entropy or SjSi, then the variations 
A/i* cannot be arbitrarily large. This is also the case when one has convergent Gaussian 
integrals. However one typically has highly peaked distributions, and thus any physical 
cutoff is naturally larger than the width of the Gaussian. Therefore a good approximation 
which yields a simple analytic result is given by integrating over the whole real line, which 
is what we always do. In our case, the variations in AM and AJ have to be smaller 
than a certain quantity Am,j such that: a) they are much smaller than the total mass and 
angular momentum, b) they are small enough to be considered independent (since for large 
values they would be related the extremality bound x < 1), and c) the truncated series 
expansions ()4.7|) and ()4.8|) remain good approximations. The latter condition can be used 
to give a simple estimation of the cutoff near extremality, namely: 

Am,j ~ 0{e) , (A.2) 

for e = 1—x. This is obtained by computing the third derivatives of the equilibrium entropy, 
which behave as ~ e~5/2 near the vertical asymptote. These give a simple estimation of 
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the order of magnitude of the next contribution in the expansion of both S and 6(3i. 
Making them neghgible w.r.t. the second order contribution (see Eqs. fj3.18p and ()3.20|) ) 
yields fTOIl . 

In our case we always have one positive and one negative eigenvalue. Now, depending 
on which quantity we want to compute, we can always arrange the integrand in such a way 
that the contribution arising from the positive ( "divergent" ) mode cancels out in the final 
result for any choice of the cutoff Am,j- As a matter of fact, one can check that: 

{Am = 00, Aj) = -Hmm + P\ 

(A.3) 

{5uj5uj) {Am, Aj = 00) = -Hjj + u^. 

This results are exact and do not require any explicit form of the coefficients in the ex- 
pansion of the entropy, only that Hmm < 0, Hjj < and detH < 0, as in our case. 
The final dependence in the cutoff that we retain finite in each case (which would give a 
divergence if taken to be infinite) exactly cancels in the numerator and the denominator 
of the general expression (jA.lll . To get this result we are approximating, in each case, the 
integral over the negative ("convergent") mode with a finite cutoff with the integral over 
the whole real line. To justify such an approximation, let us note that the width cTj of the 
Gaussian integrals is given by {i = M, J) 

ai^l/^^i^O{€"^) (A.4) 

near extremality (see Eqs. ()3.18|1 and ()3.2()j) ). This is always smaller than the cutoff ()A.2j) . 

Finally, the easiest way to find a simple approximation for the mixed correlations 
{5(3 5uj) is the following. One can check that: 

{{(3 + 5^) 5^;) (Am = 00, Aj) = -Hmj , (A.5) 

while, on the other hand: 

{6tu) {Am, Aj = 00) = . (A.6) 
Since {{(3 + Sp) SQ) = j3{5uj) + (5/5 5tD), one gets: 

{5d5u) = -HMj + l3uj. (A.7) 

Equations ()A.3|) and ()A.7|1 complete therefore the result ()4.9|) anticipated in Section HJ 

Of course, the calculation above can also be carried out by using the symmetric expres- 
sions: 

{{u + 5u) 6P) {Am, Aj = 00) = -Hmj , {SP) {Am = 00, Aj) = -f3 . (A.8) 

Note the shifts in the mean fiuctuations, {6(3) = —(3 and {6lo) = —00, which are due to the 
nonvanishing linear terms in the expansion of the entropy. 
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